Yielding of binary colloidal glasses 



T. Sentjabrskaja," E. Babaliari,* J. Hendricks," M. Laurati,*" G. Petekidis* and S.U. Egelhaaf a 

Received Xth XXXXXXXXXX 20XX, Accepted Xth XXXXXXXXX 20XX 
First published on the web Xth XXXXXXXXXX 200X 
DOI: 10.1039/b000000x 



o 

(N 



The rheological response, in particular the non-linear response, to oscillatory shear is experimentally investigated in colloidal 
glasses. The glasses are highly concentrated binary hard-sphere mixtures with relatively large size disparities. For a size ratio of 
0.2, a strong reduction of the normalized elastic moduli, the yield strain and stress and, for some samples, even melting of the 
glass to a fluid is observed upon addition of the second species. This is attributed to the more efficient packing, as indicated by 
the shift of random close packing to larger total volume fractions. This leads to an increase in free volume which favours cage 
deformations and hence a loosening of the cage. Cage deformations are also favoured by the structural heterogeneity introduced 
by the second species. For a limited parameter range, we furthermore found indications of two-step yielding, as has been reported 
previously for attractive glasses. In samples containing spheres with more comparable sizes, namely a size ratio of 0.38, the cage 
seems less distorted and structural heterogeneities on larger length scales seem to become important. The limited structural 
changes are reflected in only a small reduction of the moduli, yield strain and stress. 
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1 Introduction 

Many particle dispersions used in applications, for example 
paint, ink, cement, ceramics or foodstuffs, are characterised by 
a size distribution of the dispersed phase. Even if a monodis- 
perse system is desirable, it is often difficult to avoid a distri- 
bution of particle sizes. Furthermore, through the size distri- 
bution, the properties of a dispersion, such as its rheological 
behaviour, can be tuned, for instance to meet processing or 
application needs. To investigate the effect of a distribution of 
sizes, binary mixtures of spherical colloidal particles represent 
the simplest model system. 

The interactions and the phase behaviour of binary colloidal 
hard-sphere mixtures have been studied by theoryW£land sim- 
ulations.EKl j n equilibrium, binary colloidal mixtures exhibit 
a wider fluid-solid coexistence region than one-component 
systems, which has been thoroughly investigated in exper- 
iments. ^£23 Additionally, formation of complex crystalline 
structures through co-crystallisation of the two species is pre- 
dicted and observed.Q-QlHl For size ratios 8 = Rs/Rl % 0.2, 
where R$ and Ri are the radii of the small and large spheres, 
respectively, theory expects fluid-fluid and solid-solid coexis- 
tences, SI which are also observed in simulations^ but not yet 
in experiments. In addition, non-equilibrium glass states have 
been predicted theoreticallyL^U^] and observed experimen- 
tallyJ2122lln particular, Mode Coupling Theory (MCT) predicts 
that, at constant total volume fraction 0, a one-component 
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glass is melted upon addition of a sufficient amount of spheres 
with a different size (<5 < 0.65). 16 17 This is consistent with the 
faster structural relaxation experimentally observed in sam- 
ples with 5 « 0.6, (j) « 0.58 and intermediate mixing ratios.^ 
This leads to a strong decrease of the viscosity, which has been 
determined in experiments and simulations for a sufficiently 
large degree of mixing.'^^l Recent MCT results 17 further- 
more predict that for a large size disparity, 5 < 0.2, differ- 
ent glass states exist, which are distinguished by caging of 
one or both species, or by depletion induced bonding of the 
large spheres. The latter, for which some experimental ev- 
idence exists for 8 s» 0.1, 9 is expected to show similarities 
with attractive glasses as those observed in colloid-polymer 
mixtures.! 23 ! 24 ! 

Similar to the interactions and the phase behavior, also the 
rheological response of binary mixtures changes upon vary- 
ing the si ze an d mixing ratios. This has been studied exper- 
imentallyEiEIK] theoretically'^EI and by simulations.^^!] 

In the granular limit, i.e. when Brownian motion becomes 
irrelevant, binary mixtures with a size ratio 5 = 0.2 exhibit 
a minimum of the viscosity at a relative volume fraction of 
small spheres, xs ~ 0.4, which is known as Farris effect. 
In contrast, for colloidal mixtures a minimum of the viscos- 
ity is only observed at high total volume fractions (j) > 0.4 
and at a mixing ratio which depends on and 5. ^ With de- 
creasing 8, the minimum occurs at smaller fractions of small 
spheres, which results from a balance between the more effi- 
cient packing, since small spheres can fill the space between 
large spheres, and the depletion attraction induced between 
large spheres.' 3 - 4 -' Nevertheless, the rheology of concentrated 
binary colloidal mixtures has hardly been studiedJ^U 3 - 9 -! espe- 
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cially of spheres with significantly different sizes, i.e. small 
size ratios 8. 

Here we investigate the rheology of dispersions contain- 
ing binary mixtures with small size ratios, 8 s» 0.2 and 0.38, 
over a broad range of total volume fractions and mixing 
ratios, characterized by the relative volume fraction of small 
spheres xs = (fe/Q. Their response to oscillatory shear is 
studied with a particular focus on the non-linear viscoelas- 
tic properties, while the linear response, together with the 
structure and d ynamics at rest, will be discussed in detail 
elsewhere j 4Q l 41 l In the present case of spheres with signifi- 
cantly different sizes (i.e. small 8), the non-linear response 
contains contributions related to the different length scales 
present in the samples. This is similar to colloid-polymer mix- 
tures, where systems with attractive interactions, such as gels 
or attractive glasses, are characterized by two yielding pro- 
cesses. 42 45 The two yielding processes reflect the breaking 
of inter-particle 'bonds' and cluster breaking, in the case of 
gels, or irr eversi ble cage deformation, in the case of attrac- 
tive glasses The yielding behaviour of attractive systems 
is hence different from the one of repulsive systems, which 
typically only show one yielding mechanism related to cage 
distortion.lSlSlSE3 

2 Materials and Methods 
2.1 Rheology 

Rheological measurements were performed with an 
AR2000ex stress-controlled rheometer, and ARES G2 
and ARES strain-controlled rheometers from TA instruments, 
using cone and plate geometries of diameter D = 20 mm, 
cone angle a = 2° and gap d = 0.054 mm (AR2000ex), 
D = 25 mm, a = 2° and d = 0.048 mm (ARES G2) and 
D = 25 mm and 50 mm, a = 2° and d = 0.048 mm (ARES). 
Solvent traps were used in all rheometers to minimize 
solvent evaporation. The temperature was set to T = 20 °C 
and controlled within ±0.1 °C via a standard Peltier plate 
(AR2000ex, ARES) or an advanced Peltier system (ARES 
G2). The effects of sample loading and aging were reduced 
by performing the following rejuvenation procedure before 
each test. Directly after loading, a dynamic strain sweep 
was performed to estimate the strain amplitude y at which 
the system starts to flow, i.e. oscillatory shear was applied 
to the samples with frequency CO — 1 rad/s and increasing 
y until the sample showed a liquid-like response. Then, 
before each measurement, flow of the sample was induced by 
applying oscillatory shear at a sufficiently large strain. In the 
case of the size ratio 8 — 0.20, y = 300% was used for all 
samples. For 8 = 0.38, different values 200% < y < 1000% 
were used depending on the volume fraction and relative 
volume fraction of small particles x s . Shear was applied 



until a steady-state response, i.e. a time-independent storage 
G' and loss modulus G", was achieved, which typically 
took about 200 s. Subsequently the samples were sheared 
at 0.1% < y < 1.5% (depending on sample) until the lin- 
ear viscoelastic moduli reached a time-independent value, 
typically after 100 s to 900 s (depending on sample). This 
indicated that no further structural changes occurred and 
hence a reproducible state of the sample was reached and a 
new measurement could be started. Note that ageing effects 
might be present at longer waiting times. Measurements 
with serrated and smooth geometries, respectively, yielded 
comparable results suggesting the absence of wall slip. 

2.2 Samples 

Polymethylmethacrylate (PMMA) spheres sterically stabi- 
lized with a layer of polyhydroxystearic acid (PHSA) were 
dispersed in a mixture of cycloheptyl bromide (CHB) and 
cis-decalin that closely matched the density and refractive 
index of the colloids (8 = 0.20 and 0.19) or in a mixture 
of octadecene and bromonaphtalene which minimizes sol- 
vent evaporation (8 — 0.38). For samples in octadecene- 
bromonaphtalene, measurements of the time evolution of 
the linear viscoelastic moduli indicate the absence of sig- 
nificant gravitational effects over times much longer than 
typical measurement times. In the CHB/decalin mixture, 
the particles acquire a small charge which was screened by 
adding 4 mM tetrabutylammoniumchloride. 48 In this case, 
the colloids behave hard-sphere-like in both solvent mixtures. 
PMMA spheres with different average radii were used; R? = 
880 nm (polydispersity 0.057) and R^i = 175 nm (polydis- 
persity 0.150) to result in 8 = 0.20, R^ F = 942 nm (polydis- 
persity 0.06) and the same /?si to result in 8 = 0.19, /?l = 
358 nm (polydispersity 0.140) and Rsi = 137 nm (polydis- 
persity 0.120) to result in 8 — 0.38. The radii and polydis- 
persities were determined from the angular dependencies of 
the scattered intensity and the diffusion coefficients, obtained 
using static and dynamic light scattering, respectively, with 
very dilute colloidal suspensions (0 ~ 10~ 4 ). For the large 
spheres, a similar radius, /?£ = 885 nm, has been estimated 
from the position of the first peak of the radial distribution 
function, which was obtained by confocal microscopy.! 4 -^ Con- 
focal microscopy could be performed with these large spheres, 
because they were fluorescently labelled with nitrobenzoxa- 
diazole (NBD). Confocal microscopy was also used to de- 
termine the volume fraction of a dispersion of these spheres 
as follows. A random close packed sample was obtained by 
sedimenting a dilute suspension in a centrifuge. The sed- 
iment, whose volume fraction was roughly estimated using 
simulation results,^ was subsequently diluted to a volume 
fraction ~ 0.4 and imaged by confocal microscopy. The 
imaged volume was partitioned into Voronoi cells and their 



2 



10 



mean volume determined. The ratio of the particle volume 
to the mean Voronoi volume provides an estimate of the vol- 
ume fraction of the sample, = 0.43. This allowed us to cal- 
culate the volume fraction of the random close packed stock 
solution 0r CP =0.68. The smaller spheres were too small to 
be imaged (thus also not fiuorescently labelled). The volume 
fraction of their sediment was estimated taking into account 
their polydispersity:^! <^I p ~ 0.68 for spheres with radius 
Rsi = 175 nm and — 0rcp — 0-^7 for spheres with radii 
Rsi = 137 nm and /?l = 358 nm. The value of the volume 
fraction is known to suffer from relatively large uncertain- 
ties.^ Thus the value of obtained for the large spheres was 
used as a reference value and the volume fraction of the two 
batches containing the smaller particles adjusted using Theo- 
logical measurements as follows. Linear viscoelastic moduli 
for samples at nominally equal volume fraction (0 ~ 0.58 for 
8 = 0.20, ~ 0.595 for 8 = 0.38, ~ 0.61 for 8 = 0.19) were 
measured in Dynamic Frequency Sweeps (DFS) at a strain 
amplitude 0.1% < y < 1.5% (depending on sample). The ob- 
tained storage moduli G' and loss moduli G" as a function of 
oscillation frequency co are expected to agree for spheres of 
different size but the same volume fraction, if the moduli are 
rescaled by the energy density ~ IcqT j R 3 and the frequency by 
the Brownian time Tb = R 2 /Dq with Do = 6nr\R the Stokes- 
Einstein-Sutherland diffusion coefficient in the dilute limitf^! 
and T] the solvent viscosity. The dispersion of small spheres 
was diluted until its rescaled linear response matched that of 
the dispersion of large spheres with the desired volume frac- 
tion, i.e. until an equivalent rheological response in the linear 
regime was obtained (Fig.[TJ. Furthermore, it was verified that 
the normalised elastic modulus G' and its dependence coin- 
cides, for all particles used, with that of a dispersion contain- 
ing crystallising colloids with a low polydispersity, whose vol- 
ume fraction was determined in the crystal-fluid coexistence 
region. 52 If illuminated by laser light, Bragg reflections were 
not observed, indicating the absence of crystallinity in the one- 
component dispersions. By mixing appropriate amounts of the 
one-component dispersions, samples with different total vol- 
ume fractions and relative volume fractions of small parti- 
cles x s = 0s/ (0s + 0l) were prepared, where 0s and 0l are 
the volume fractions of small and large particles, respectively. 
Samples with different xs and two different values of (for 
8 = 0.20, 0.38) as well as fixed x$ = 0.65 and different values 
of (for 8 = 0.19) were investigated. 

3 Results and Discussion 

In Dynamic Strain Sweep (DSS) experiments, a sinusoidal 
strain is applied whose frequency CO is constant but whose 
amplitude y is increased in steps, starting in the linear vis- 
coelastic regime and progressing into the non-linear regime. 
The stress response of the system is recorded as a function 
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Fig. 1 (top) Storage, G' (full symbols), and loss, G" (open symbols), 
moduli of samples containing large (□) and small (O) spheres, 
respectively, as a function of frequency co obtained by Dynamic 
Frequency Sweep measurements for (left) a size ratio S = 0.20 and 
total volume fraction <j> = 0.58 and (right) 8 = 0.38 and = 0.595. 
(bottom) Same data in units proportional to the energy density, 
i.e. kftT/R 3 , and Brownian time Zg = Dq/R 2 . The strain amplitude 
was y= 0.5% for 8 = 0.20 and 7= 1.5% for 8 = 0.38. 



of strain amplitude y. Figure [2] shows results of DSS mea- 
surements for samples with size ratio 8 = 0.20, total volume 
fractions = 0.61 and 0.58, and different relative volume 
fractions of small spheres xs- Beyond the linear viscoelastic 
regime the stress response in DSS experiments significantly 
deviates from a simple sinusoidal form and can be decom- 
posed into higher order (odd) harmonics, as shown before for 
one-component hard-sphere glasses. 52 . However, the G' and 
G" values shown in figure [2] correspond to the first harmonic 
contribution of the stress response. To allow for a compar- 
ison of the different samples, measurements were not per- 
formed at a constant frequency CO, but at a fixed oscillatory 
Peclet number Pe m — 0)Tb. It is the ratio of the Brownian 
time of the system, Xg = (R 2 /Dq), and the timescale imposed 
by shear, i.e. the inverse of the frequency, T ffl = 1/(0. Thus, 



Pe a = ((67rT7/? 3 )/(Ar B r))co and 



(R 3 )=4 



(1) 



We applied Pe m = 5.55 x 10 _1 corresponding to 7.6 x 
10~ 2 rad/s < CO < 9.7 rad/s, depending on Xs- 

The one component systems (x$ = and 1) for both 
show th e charact eristic response of a hard sphere glass 
(Fig. [2^,b)J^ & 3 l 47 l 52 l (Note that due to the much lower energy 
density of the samples with the large spheres, their response 
is much weaker and thus more affected by noise.) The stor- 
age modulus G' is larger than the loss modulus G" in the lin- 
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Fig. 2 (a,b) Storage, G', (full symbols) and loss, G", (open symbols) moduli as a function of the strain amplitude y obtained in DSS 
measurements. The size ratio 8 = 0.20, the total volume fraction (a) <f> = 0.61 and (b) <j) = 0.58, the relative volume fraction of small particles 
X S = 0.0 (i ), 0.1 (A), 0.3 (V), 0.5 (0), 0.7 ( l), 0.9 (O), 1.0 (□) and Pe m = 5.55 x 10" 1 (corresponding to 7.6 x 10~ 2 rad/s < a < 9.7 rad/s). 
Circles indicate the yields strains y y and the red solid line their jrg -dependence. 



ear viscoelastic regime, with their values comparable to the 
ones obtained in dynamic frequency sweeps (Fig. [T|l.E2l The 
two moduli become equal at a strain amplitude y y (highlighted 
with circles in Pig. |2j, which is identified with the yield strain 
of the glass. At the yield strain y y and the corresponding yield 
stress Oy, the local environm ent of a pa rticle is irreversibly re- 
arranged, i.e. its cage broken .I42J43J46J47J p or y > y qii j s j^-gg]- 

than G' and the system starts to flow. In this regime, G" shows 
a maximum which indicates the largest energy dissipation and 
has previously also been used to estimate the yield strain asso- 
ciated with irreversible rearrangements of the cage 

PHUpon 

increasing the volume fraction from = 0.58 to 0.61, the lin- 
ear viscoelastic moduli and the yie ld strain y y increase. This 
is consistent with previous studies, ^ 47 ! 52 ! which found y y to 
increase with volume fraction up to « 0.62, beyond which it 
decreases due to the approach toward random close packing. 

Keeping the total volume fraction (j) constant, but chang- 
ing the composition to xs = 0.9, the storage and loss moduli 
decrease (Fig. [2j. The decrease is not only due to the pres- 
ence of large particles and hence a lower energy density, but 
remains even if the moduli are rescaled by the energy density 
(nkftT) ~ l/(7? 3 ). This indicates a softening of the glass. A 
softer response is also reflected in a reduced yield strain y y 
and yield stress rj y (circles in Fig. [2] A further decrease of the 
relative volume fraction of small spheres to xs = 0.7 leads to 
an additional reduction of the storage, G', and loss, G", modu- 
lus, yield strain y y and stress rj y , which indicates that the glass 
still becomes mechanically weaker. Forxs > 0.7, comparable 



effects are found for <j) = 0.61 and 0.58. 

This is different for xs < 0.7. For the higher total volume 
fraction = 0.61, the samples with xs = 0.5 and 0.3 have a 
much smaller G' which, however, is still slightly larger than 
G" and the samples hence show a weak solid-like response in 
the linear viscoelastic regime (Fig. [2^). This is consistent with 
Yy and rj y values which are more than one and almost three or- 
ders of magnitude smaller, respectively, than typical values of 
one-component hard-sphere glasses at the same total volume 
fraction. Hence the samples become very brittle and may flow 
plastically at smaller strain amplitudes or stresses. A closer 
inspection of the response of the sample with <p = 0.61 and 
xs = 0.3 reveals a particularly interesting strain amplitude (y) 
dependence of the moduli (Fig. [3]). The linear response ends 
already at ym 0.2% (Fig. [3] arrow on the left), beyond which 
G' decreases smoothly up to y ps 4%, where it shows a kink 
and subsequently decreases with a power-law, while G" shows 
a small maximum (Fig. [3] arrow on the right). This response 
suggests the presence of two length scales, most likely asso- 
ciated with the small and large spheres, which both contribute 
to the yielding of the system at this xs. The first yielding at 
small strains ym 0.2% might correspond to plastic rearrange- 
ments of cages formed by small spheres. Cage distortion and 
yielding might be facilitated by the shear-induced interaction 
with the large spheres, i.e. contact forces between large and 
small spheres. Once these cages are rearranged, the system 
is still prevented from flowing by the cages of large spheres 
which are only slightly deformed. At strains of about 4% the 
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Fig. 3 Storage, G', (full symbols) and loss, G", (open symbols) 
moduli as a function of the strain amplitude y obtained in DSS 
measurements for S = 0.20, a relative volume fraction of small 
spheres xs = 0.3 and the total volume fraction = 0.61 (O) and 
<j) = 0.58 (□). Arrows indicate the two yielding points observed for 
the sample with (j) = 0.61. 

cages of large spheres deform and the system starts to flow. 
The ratio between the two yield strains, 0.2/4 ~ 0.04, cor- 
responds to 8/4 which suggests a non trivial scaling of the 
yield strains with the cage size (which would give a factor 8). 
This finding could also result from the moderate polydisper- 
sity of the small spheres, which implies a distribution of the 
effective size ratio and accelerates the dynamics,^ and could 
contribute to smear out the double yielding phenomenon. A 
two-steps yielding behavior has also been observed for attrac- 
tive glasses and gels.^MSICompared to (j) = 0.61, for the lower 
total volume fraction <p = 0.58 decreasing the relative volume 
fraction of small spheres to xs = 0.5 has an even stronger ef- 
fect (Fig. |2j3). Within the whole examined range of strain am- 
plitudes, G" > G' implying fluid-like behavior. Thus, the glass 
is melted. Fluid-like behavior in the whole range of measured 
y is also observed for xs = 0.3, with the response being sim- 
ilar to that obtained for (j) =0.61, except for the smallest y 
(Fig. |3J. Samples showing fluid-like behavior (xs = 0.3, 0.5) 
do not present a finite value of yield strain and stress, corre- 
sponding to the missing circles in Fig. [2]?. The melting of the 
glass is caused by the larger free volume fraction created by 
the presence of small spheres, as will be discussed in more de- 
tail later. This is similar to the behaviour of one-component 
systems when <p is decreased below the glass transition. 

Finally, the samples at both total volume fractions show the 
response of a weak solid for xs = 0.1. For <p = 0.61 the stor- 
age modulus G' is further reduced and becomes similar to G", 
indicating the proximity of a transition to the fluid state. On 
the other hand the yield strain y y and stress a y are slightly in- 
creased (Fig. [2j. In contrast, for (j) = 0.58 the response again 
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Fig. 4 (a,b) Storage, G', (full symbols) and loss, G", (open 
symbols) moduli as a function of strain amplitude y obtained in DSS 
measurements. The size ratio 5 = 0.38, (a) the total volume fraction 
<p — 0.615 and the relative volume fraction of small particles 
x s = 0.0 (>), 0.08 (A), 0.25 (0), 0.5 ( ), 1.0 (□) and (b) (j) = 0.595 
andx s = 0.0 (>), 0.25 (0), 0.5 ( l), 0.75 (O) and 1.0 (□). The 
frequency a = 1 rad/s. Circles indicate the yield strains y y and the 
red solid line their xs -dependence. 



changes qualitatively, which implies a reentrant behavior; the 
melting and re-formation of a solid glass state as the fraction 
of small spheres is reduced. 

A second size ratio, 8 = 0.38, was investigated also at two 
total volume fractions = 0.595 and 0.615 and different rel- 
ative volume fractions of small particles xs (Fig. |4j. Starting 
from the one-component systems and increasing the amount 
of the second component, the storage modulus G' decreases 
in the linear viscoelastic regime indicating a softening of the 
glass, similar to the findings with 8 = 0.20 (Fig. |7^). How- 
ever, in the case of 8 = 0.38, the minimum of G' is located at 
xs ~ 0.5 for both (j). Note that in terms of the relative number 
of small spheres ^$ = nz,/(ns + = *s[<5 3 +x^(l — <5 3 )]~\ 
where «s and «l are the number densities of small and large 
spheres, respectively, the minimum of G' is found for both size 
ratios at values of > 0.85. Furthermore, the minimum in the 
xs -dependence is much weaker for the yield stress a y and ab- 



1-015 



10 



10 



10 



5 = 0.19, x s =0.65. 
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Fig. 5 (a) Storage, G', (full symbols) and loss, G", (open symbols) 
moduli as a function of strain amplitude y obtained in DSS 
measurements for (a) The size ratio 5 = 0.19, the relative volume 
fraction of small particles xs = 0.65, total volume fractions 
= 0.61 (O), 0.59 (□), 0.57 (A) and 0.55 (0). (b) 8 = 0.38, 
x s = 0.5, <p = 0.615 (O) and 0.595 (□). Frequency CO = 1 rad/s, 
corresponding to Pea, = 8.99 x 10~ 2 for 8 = 0.19 and 
Pe a = 8.85 x 10~ 2 for 8 = 0.38. (c) Yield stress oy, (d) yield strain 
Y y and (e) ratio G'/G" in the linear viscoelastic regime (y = 0.5% 
and 1% for 5 = 0.2 and 0.38, respectively), as a function of total 
volume fraction <j> for samples of plot (a) (O) and (b) (□). 



sent for the yield strain y y (Fig. |7j?,c). Thus, no melting of the 
glass is observed for 8 = 0.38. 

Having studied the rheological response as a function of the 
relative volume fraction of small spheres xs, we now turn to 
the dependence on the total volume fraction (j) for constant 
x s = 0.65 (8 = 0.19) and 0.5 (8 = 0.38) (Fig. |). With de- 
creasing 0, the storage modulus G' decreases in the linear 
regime and approaches the loss modulus G" (Fig. [5^,b, e, 
which shows the ratio G'/G"). Thus, with decreasing 0, the 
solid-like response becomes weaker. This is particularly pro- 
nounced for 8 = 0.19, which shows a fluid-like response for 
<p = 0.55, that is G" > G' in the linear viscoelastic regime 
(Fig.|5]a,e). The yield point, i.e. the yield strain y y and stress 
Oy, decreases with decreasing (j) for both values of 8 and, for 
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Fig. 6 Storage, G', (full symbols) and loss, G", (open symbols) 
moduli as a function of strain amplitude y obtained in DSS 
measurements, (left) The size ratio <5 = 0.20, total volume fraction 
(top) (j) = 0.61 and (bottom) 0.58, relative volume fraction of small 
spheres xg = 0.9, and frequencies CO = 1 rad/s (0), 5 rad/s (<) and 
10 rad/s (O). (right) 8 = 0.38, (top) <j> = 0.615 and (bottom) 0.595, 
x s = 0.5, CO = 0.1 rad/s (V), 1 rad/s (0) and 10 rad/s (O). 



8 = 0.19 disappears at = 0.55, i.e. the sample becomes a 
fluid (Fig. |5};,d). This is consistent with the response of one- 
component systems, whose yield strain y y also decreas es with 
decreasing until a transition to a fluid occurs .E § l 47 l 52 l 

The decrease in y y is attributed to the fact that, upon de- 
creasing 0, the cages become larger and looser and thus in- 
creasingly smaller distortions of the cages are sufficient to 
allow the particles to escape through Brownian motion. Fi- 
nally, in the fluid phase (if) = 0.55), particles can leave the 
cage even in the absence of shear. The sample with 8 = 0.19, 
xs = 0.65 and <p = 0.57 shows a dependence of G' and G" 
on the strain amplitude y similar to that of the sample with 
8 = 0.20, xs = 0.3 and = 0.61 (Fig. [3j, which again sug- 
gests the presence of two yielding points. Note that this sam- 
ple is a dense, slowly relaxing fluid and not a glass, according 
to the frequency dependence of the linear viscoelastic moduli 
(data not shown). Nevertheless, the similarity of the response 
of the two samples suggests that a glass state similar to that of 
8 = 0.20, xs = 0.3 and (j) = 0.61, i.e. characterized by a dou- 
ble yielding process and caging on two lenght scales, might be 
obtained at slightly larger for xs = 0.65. This is in agree- 
ment with MCT predictions, where a transition from a glass 
characterized by caging on one length scale (that of the small 
spheres) at high 0, to a glass characterized by caging on two 
length scales at lower 0, and successive melting of this glass 
wi th fur ther decreasing <p, is expected along a line of constant 

x s pm 
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The results of the DSS measurements show a slight depen- 
dence on frequency (Fig.|6]l. In the linear viscoelastic regime, 
the storage modulus G' increases with increasing frequency 
CO, in agreement with results of our Dynamic Frequency 
Sweep (DFS) measurements (Fig. [T} and as discussed in 
more detail elsewhere.* With increasing frequency co, the 
probed times decrease and are progressively shorter than 
the structural relaxation time. This leads to an increasingly 
more elastic response. Also the yield strain y y and stress ay 
increase with increasing Pe m (Fig.[7]3,c). Thi s is simila r to the 
behavior of one-component colloidal glasses^ 42 ! 43 ! 46 ! 47 ! and can 
be understood as follows. Shear-induced cage deformation 
facilitates the escape of particles from their cage through 
Brownian motion, which results in yielding. In oscillatory 
shear, the maximum cage deformation is achieved at the 
largest excursion. In the vicinity of this point a particle is 
most likely to escape from the cage by Brownian motion. 
With increasing frequency, the particles spend less time at 
the maximum (but more frequently) and are therefore less 
likely to escape because the escape probability depends rather 
on the balance between the residence time at the maximum 
and the Brownian time than on the attempt rateJ 44 l 55 ^D. The 
reduced escape probability must be compensated by a larger 
cage deformation. Thus, with increasing frequency CO, a 
larger strain and stress will be required, and hence stored, 
before the cage breaks. 

Our findings are summarized in Fig. [7] For a given to- 
tal volume fraction 0, adding a second component to the 
one-component systems results in a weaker elastic response. 
For 8 = 0.20, the glass softens particularly strongly and, if 
the sample is sufficiently close to the glass transition (here 
= 0.58), even melts, that is shows a fluid-like response. 
This reduction in G' is not symmetric with respect to the 
one-component systems, but is more pronounced for glasses 
mainly consisting of large spheres to which a small amount 
of small spheres has been added. This is evident when com- 
paring, for example, G' for samples with xs — 0.1 and 0.9. 
This asymmetry might, however, be due to the choice of the 
control parameter, here the relative volume fraction of small 
particles x$. Instead, one could use the relative number of 
small spheres, Hence xs =0.1 corresponds to = 0.93 
while xs = 0.9 implies a relative number of large spheres of 
only = 8.9 x 10~ 4 . This might explain why for xs = 0.9 
the cage of small spheres is not significantly affected by the 
small number of large spheres. In contrast, for x§ = 0.1 a 
large number of small spheres has to be accommodated by 
the large spheres, which is likely to induce a significant cage 
deformation and to result in a significant softening. This is 
supported by confocal microscopy measurements of the struc- 
ture and dynamics of the large spheresJ^ED which are sum- 
marized in Figj7jl,e. Already at x$ = 0.1, the pair distribu- 
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Fig. 7 (a) Storage modulus G' in the linear viscoelastic regime 
(7 = 0.2%), (b) yield strain y y estimated from the crossing point of 
G' and G" and (c) corresponding yield stress o~ y as a function of the 
relative volume fraction of small particles xs for samples with size 
ratio 8 = 0.20 and total volume fraction (j) = 0.61 (■) and <j> — 0.58 
(•) and Pea, = 5.55 10" 1 , and 8 = 0.38, (?) = 0.615 (□) and 
(j) = 0.595 (O) and CO = 1 rad/s. (b) also contains results for 
Pe m = 2.75 10" 1 (x) and Pe a = 5.55 10~ 2 (+) for the sample with 
xs = 0.9 and <j> = 0.61. (d) Height of the maxima of the pair 
distribution function g(r), g ma x, corresponding to r = 2Ri (o), 
r = 2(R L + R S ) (♦) and r = 2(R L + 2R S ) (*) and (e) localisation 
length L extracted from plateaus of mean squared displacements as a 
function of xs, for samples with 8 = 0.2 and <p = 0.61 Error 
bars are smaller than the symbols in all plots. 



tion function g(r) does not only show a peak at r = 2R^, but 
also a shoulder at r = 2(R\ j + Rs) indicating that the cage of 
large spheres is deformed and that a significant fraction of 
large particles is separated by small particles. This cage de- 
formation leads to a slight increase in the particle localisation 
length extracted from the plateau of mean-squared displace- 
ments, but the dynamics of the system is still arrestedi^l 4 -!! 
The reduced localisation is thought to be responsible for the 
strong decrease in yield strain. Rearrangement of the cage 
of large spheres becomes even more pronounced as xs is in- 
creased to 0.3 and 0.5, as demonstrated by the increasingly 
larger reduction of the peak at r = 2/?l an d the corresponding 
increase at r = 2 (i?L + ) a s well as the appearance of ad- 
ditional peaks at distances r = 2(R^ + nRs,), with n an integer 
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number. For these xs, the dynamics show diffusion at <j) = 0.58 
and sub-diffusion at <p = 0.61 with a decreasing localization 
length suggesting that large particles start to be localized more 
tightly by small spheres.^HUDuring this process of cage rear- 
rangement, for = 0.58 the elasticity decreases and the yield 
strain y v disappears due to the melting of the glass, while at 
= 0.61 both G' and y y start to increase again significantly 
above x$ — 0.3 possibly due to the emergent caging and lo- 
calisation of large spheres by small spheres. For x$ > 0.5 the 
localisation in cages of small spheres, i.e. the transition to a 
different glass state, is accomplished: Large particles repre- 
sent a dilute phase in a dense matrix of small spheres and are 
localised on distances which are about a factor 8 = 0.2 smaller 
than at x$ — and their dynamics are again arrested j 4Q l 41 l The 
tighter localisation and dynamical arrest induce an increased 
G' and y y towards the values of the one-component glass of 
small spheres. A pronounced effect of size and mixing ratios 
on the structure and dynamics of the glass was also reported 
for 2D colloidal glass formersJSKI In particular, changes in 
the relative content of the small component and the size ratio 
ha ve bee n reported to have pronounced effects on the dynam- 

For 8 = 0.38, the softening is less pronounced and no melt- 
ing is observed. Moreover, the dependence of G' on x§ is more 
symmetrical with respect to the one-component systems. The 
smaller effect is attributed to the fact that the small particles 
have a reduced ability to occupy the interstitial space between 
the large particles at this size ratio. The critical value 8 C at 
which the small spheres cannot fill the space in between two 
large ones in a dense packing of large spheres can be esti- 
mated: In a group of 9 spheres arranged as in a body-centred 
cubic lattice and in contact with each other, the centers of two 
spheres along a face diagonal are separated by 2\f2R\ j and a 
small sphere can fill the space left in between the large spheres 
if R s < (V2- l)R L w 0.41fl L , i.e. 8 < 8 C » 0.41, which is 
comparable to 8 = 0.38. Although in the glass states consid- 
ered here, ordered configurations are not expected, the size of 
the void space might be similar. Thus the cage itself, i.e. the 
first neighbo r shell, is not expected to be rearranged signifi- 
cantly and the softening hence appears to be caused by the het- 
erogeneity of the cage on an intermediate length scale rather 
than a more efficient packing. The weaker cage deformation 
induced by the smaller packing ability at 8 — 0.38 can also 
explain the weaker reduction of the yield strain and stress ob- 
served at intermediate mixing ratios for this 8. 

Instead of the relative volume, Xs, or number, 4s, fraction 
of small particles, we now consider the distance to the two 
limiting volume fractions of the glass state, corresponding 
to the glass transition and random close packing. Mode 
Coupling Theory (MCT) predicts 17 that in mixtures the glass 
transition is shifted to higher total volume fractions. For 
example, for the size ratio 8 — 0.38 the maximum volume 
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Fig. 8 (a) Storage modulus G' in the linear viscoelastic regime and 
(b) yield strain y y as a function of the free volume 0f ree for the same 
samples as presented in Fig. [7] and for samples of Fig.[5]( ). The 
red line in (a) shows a power-law fit G'{/? 3 ) /ksT ~ 0^3, with P~ 3. 
The dashed lines represent data of one-component hard-sphere 
glasses.^! 



fraction for the glass transition is expected at xs ~ 0.4, which 
is consistent with the occurrence of maximum softening in 
our experiments. The shift of the glass transition could be 
related to the addition of small particles with their larger 
mobility. This might favour structural rearrangements of 
the large spheres through collective motions and lead to a 
glass with a reduced elasticity, i.e. G'. In mixtures, MCT 
predicts qualitative changes of the relative particle mobilities, 
associated with different glass states. In addition, the more 
efficient packing in mixtures results in an increased total 
volume fraction at random close packing, 0rcp- Theoretical 
predictions for 0rcp are available for binary mixtures of 
monodisperse hard spheres, with different size ratios 8 and 
mixing ratios, i.e. xs. 6263 Based on the predictions for 0rcp> 
we calculate the available free volume 0f ree = 0rcp — as a 
function of xs and 8. (Predictions for 8 = 0.17 and 0.39 are 
used for the experimental 8 = 0.20 and 0.38, respectively.) 
Note that predicted values of 0rcp where shifted by the 
difference between the value of 0rcp in the monodisperse 
case (0rcp = 0.4) and the experimental values of 0rcp 
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(0rcp = 0.68 for 8 = 0.2 and R cp 
section P 



0.67 for 8 = 0.38, 



2.2 1. With decreasing free volume 0f ree , that is toward 



random close packing, the storage modulus G' is found to 
increase (Fig. IsJ. The dependence of G' on 0f ree indicates 
a common behavior for all 8 and investigated and can 
be approximately described by a power-law dependence 
G'(R 3 )/kgT ~ (0free)~ P i with p 3. A similar power-law 
dependence is observed for one-component hard- sph ere 
systems up to </> free < 0.1 (Fig. § dashed line).' 46 ' 47 ' 52 ' ^ At 
larger values of </>f ree the one-component system shows a 
sharper decay. 



We now consider the dependence of the yield point on the 
free volume 0f ree . The x$ -dependence of the yield strain y y and 
stress <7 y is quite different for the two size ratios (Fig.|7j3,c). In 
particular, both, y y and cr y , show a much weaker dependence 
on xs for 8 = 0.38 than 8 = 0.2. This can also be linked to 
the free volume available for structural rearrangements. The 
dependence of y y on 0f lee (Fig.|8j3) indicates that toward small 
free volumes, the yield strain saturates at an approximately 
constant value y y « 20%, which ag rees with the yield strain 
observed in one-component glasses ffi l 43 ! 52 ! 64 ! At smaller val- 
ues of 0f ree , i e. very close to RCP, which are not reached 
here, in the one-component systems the yield strain decreases 
(Pig j8]>. In contrast, toward large 0f ree > 0.1 a strong decrease 
of Yy is observed (for samples with 8 = 0.20 since only they 
reach large enough 0f ree due to their large 0rcp). This de- 
crease indicates that if a sufficiently large free volume, i.e. a 
sufficiently loose packing, is present, significant structural re- 
arrangements can be induced by small strains. Their strong de- 
crease of the yield strain is observed for samples in which the 
small spheres occupy the free space between the large spheres. 
The intercalation of small spheres in between large spheres 
possibly induces a strong deformation of the cage. Similar 
effects have been observed in mixtures of star polymers with 
significant size disparity. 65 -l This supports our previous find- 
ing that yielding is not only facilitated by the increase of free 
volume but also by structural heterogeneities leading to cage 
deformation. Interestingly, the strong decrease in the yield 
strai n y y for ( fa ree > 0. 1 is not observed in one-component sys- 
t ems J46|47l52] smce m m j s regime the glass is melted. This is 
also consistent with G' sharply decreasing for 0f ree > 0.1 for 
the one-component system (Fig. [8^, dashed line). We specu- 
late that in the glass state G'(R 3 ) /k B T — (0f re e)~ p with p « 3 
for one and two-component systems. The slight shift between 
our system (red line) and the previous one-component data 
(dashed line) is due to different interactions mediated by dif- 
ferent solvents.^ These findings show that at large values 
of 0free, a glass can still be formed in the mixture (possibly 
due to attractions) while a dense fluid is observed in the one- 
component system. 



4 Conclusions 



The linear and non-linear response to oscillatory shear has 
been studied in concentrated binary hard-sphere mixtures with 
large size disparities, 8 « 0.20 and 0.38. In the linear regime, 
the response of mixtures is softer than that of the correspond- 
ing one-component systems at the same total volume fraction 
(j), as demonstrated by the smaller normalised storage modulus 
G' . The softening is associated with a shift of random close 
packing to larger total volume fractions, and thus a larger free 
volume fraction 0f ree , which results from the more efficient 
packing in two-component systems. 62 Pronounced softening 
occurs for the size ratio 8 = 0.20 and for samples containing 
a majority of large spheres (x$ < 0.5). This indicates that soft- 
ening is not only a result of an increased free volume 0f ree but 
also of cage distortions due to small particles filling the space 
between the large spheres. In contrast, in the samples with 
a smaller size disparity (8 = 0.38) and a majority of small 
spheres (x§ > 0.5), we can speculate that on average the cage 
structure should be poorly affected due to the reduced abil- 
ity of the small component to fill space in between the large 
spheres, and heterogeneities are thus only introduced beyond 
the first neighbor shell, which results in a weaker softening of 
the glass. 

In the non-linear regime, the more efficient packing in the 
mixtures affects the yielding behaviour. If the free volume 
0f ree is only slightly increased, yielding is characterised by a 
one-step cage break-up, as in one -component systems. With 
increasing free volume, yielding occurs at smaller deforma- 
tions. Interestingly, at large values of the free volume, the 
presence of a small but finite yield strain indicates the persis- 
tence of a weak solid-like state in the mixtures, while at com- 
parable free volume a one-component system melts. This oc- 
curs in systems where the small spheres can occupy the space 
in between the large spheres, which suggests that the intercala- 
tion of small spheres induces a strong deformation and loosen- 
ing of the cage structure and thus contributes to the reduction 
of the yield strain. Moreover, the yielding behaviour could be 
affected by a possible transition between different glass states, 
in particular if it is associated with the mobility of the small 
spheres, which could facilitate yielding through collective mo- 
tions. In addition to the one-step yielding behaviour, we also 
found indications of a more complex two-step yielding behav- 
ior for a sample withxs =0.3, <p =0.61. The two steps could 
be linked to the two different length scales present in these 
samples, representing caging of small and large spheres, re- 
spectively. While two length scales are present in all mixtures, 
in most samples one of the two dominates, rendering the sec- 
ond yielding insignificant. 



1-019 



5 Acknowledgments 

We acknowledge support from the Deutsche Forschungsge- 
meinschaft through the FORI 394 Research unit and EU fund- 
ing through the FP7-Infrastructures 'ESMI' (CP&CSA-2010- 
262348). We also thank A. B. Schofield for the synthesis of 
the PMMA particles and N. Koumakis for his help at the ini- 
tial stages of experiments at FORTH. We thank in addition Th. 
Voigtmann, K. J. Mutch, P. Chauduri, J. Horbach, M. Fuchs, 
R. Castaneda Priego and W. C. K. Poon for stimulating dis- 
cussions. 

References 

1 S. Amokrane, A. Ayadim and J. Malherbe, J. Chem. Phys., 2005, 123, 
174508. 

2 R. Roth, R. Evans and S. Dietrich, Phys. Rev. E, 2000, 62, 5360. 

3 D. J. Ashton, N. B. Wilding, R. Roth and R. Evans, Phys. Rev. £,2011, 
84,061136. 

4 M. Dijkstra, R. van Roij and R. Evans, Phys. Rev. E, 1999, 59, 5744. 

5 B. Gotzelmann, R. Roth, S. Dietrich, M. Dijkstra and R. Evans, Europhys. 
Lett., 1999,47,398. 

6 J. G. Malherbe and W. Krauth, Mol. Phys., 2007, 105, 2393. 

7 J. S. van Duijneveldt, A. W. Heinen and H. N. W. Lekkerkerker, Europhys. 
Lett., 1993,21, 369. 

8 A. D. Dinsmore, A. G. Yodh and D. J. Pine, Phys. Rev. E, 1995, 52, 4045. 

9 A. Imhof and J. K. G. Dhont, Phys. Rev. Lett., 1995, 75, 1662. 

10 N. Hunt, R. Jardine and P. Bartlett, Phys. Rev. E, 2000, 62, 900. 

11 P. Bartlett, R. H. Ottewill and P. N. Pusey, Phys. Rev. Lett, 1992, 68, 
3801. 

12 X. Cottin and P. A. Monson, /. Chem. Phys., 1995, 102, 3354. 

13 A. B. Schofield, Phys. Rev. E, 2001, 64, 051403. 

14 A.-P. Hynninen, L. Filion and M. Dijkstra, J. Chem. Phys., 2009, 131, 
064902. 

15 M. Dijkstra, R. van Roij and R. Evans, Phys. Rev. Lett., 1998, 81, 2268. 

16 W. Gotze and Th. Voigtmann, Phys. Rev. E, 2003, 67, 021502. 

17 Th. Voigtmann, Eur. Phys. Lett., 201 1, 96, 36006. 

18 Ph. Germain and S. Amokrane, Phys. Rev. Lett., 2009, 102, 058301. 

19 F. Tchangnwa Nya, A. Ayadim, Ph. Germain and S. Amokrane, /. Phys.: 
Condens. Matter, 2012, 24, 325106. 

20 S. R. Williams and W. van Megen, Phys. Rev. E, 2001, 64, 041502. 

21 B. E. Rodriguez, E. W. Kaler and M. S. Wolfe, Langmuir, 1992, 8, 2382. 

22 G. Foffi, W. Gotze, F. Sciortino, P. Tartaglia and Th. Voigtmann, Phys. 
Rev. Lett., 2003, 91, 085701. 

23 K. N. Pham, A. M. Puertas, J. Bergenholtz, S. U. Egelhaaf, A. Moussaid, 
P. N. Pusey, A. B. Schofield, M. Cates, M. Fuchs and W. C. K. Poon, 
Science, 2002, 296, 104-106. 

24 K. N. Pham, S. U. Egelhaaf, P. N. Pusey and W. C. K. Poon, Phys. Rev. E, 
2004, 69,011503. 

25 A. T. J. M. Woutersen and C. G. de Kruif, /. RheoL, 1993, 37, 681. 

26 P. D'Haene and J. Mewis, Rheol. Acta, 1994, 33, 165. 

27 W. Richtering and H. Muller, Langmuir, 1995, 11, 3699. 

28 W. J. Hunt and C. F. Zukoski, Langmuir, 1996, 12, 6257. 

29 P. Gondret and L. Petit, Langmuir, 1997, 41, 1261. 

30 R. Greenwood, P. F. Luckham and T. Gregory, /, Colloid Interface Sci., 
1997, 191, 11. 

31 T. Shikata, H. Niwa and Y. Morishima, J. Rheol., 1998, 42, 765. 

32 T. Ohtsuki, PhysikaA, 1983, 122, 212. 

33 G. Nagele and J. Bergenholtz, J. Chem. Phys., 1998, 108, 9893. 



34 R. A. Lionberger, Phys. Rev. E, 2002, 65, 061408. 

35 C. Y. Chang and R. L. Powell, J. Fluid. Mech, 1993, 253, 1. 

36 C. Y. Chang and R. L. Powell, J. Rheol, 1994, 38, 85. 

37 C. Y. Chang and R. L. Powell, Phys. Fluids, 1994, 6, 1628. 

38 R. J. Farris, Trans. Soc. Rheol, 1968, 12, 281. 

39 R. Foudazi, I. Masalova and A. Y. Malkin, /. Rheol, 2012, 56, 1299. 

40 T. Sentjabrskaja, M. Laurati, S. U. Egelhaaf and Th. Voigtmann, in prepa- 
ration. 

41 T. Sentjabrskaja, S. U. Egelhaaf and M. Laurati, AIP Conf. Proc, 2013, 
1518, 206. 

42 K. N. Pham, G. Petekidis, D. Vlassopoulos, S. U. Egelhaaf, P. N. Pusey 
and W. C. K. Poon, Europhys. Lett, 2006, 75, 624-630. 

43 K. N. Pham, G. Petekidis, D. Vlassopoulos, S. U. Egelhaaf, W. C. K. Poon 
and P. N. Pusey, /. Rheol, 2008, 52, 649-676. 

44 M. Laurati, S. U. Egelhaaf and G. Petekidis, / Rheol, 201 1, 55, 673-706. 

45 N. Koumakis and G. Petekidis, Soft Matter, 2011, 7, 2456-2470. 

46 G. Petekidis, A. Moussaid and P. N. Pusey, Phys. Rev. E, 2002, 66, 
051402. 

47 G. Petekidis, D. Vlassopoulos and P. N. Pusey, Faraday Discuss, 2003, 
123, 287. 

48 A. Yethiraj and A. van Blaaderen, Nature, 2003, 421, 513-517. 

49 M. C. Jenkins and S. U. Egelhaaf, Adv. Coll. Interface Set, 2008, 136, 
65-92. 

50 W. Schaertl and H. Silescu, /. Stat. Phys., 1994, 77, 1007-1025. 

51 W. C. K. Poon, E. R. Weeks and C. P. Royall, Soft Matter, 2012, 8, 21-30. 

52 N. Koumakis, A. Pamvouxoglou, A. S. Poulos and G. Petekidis, Soft Mat- 
ter, 2012,8, 427 1^1284. 

53 H. J. Schope, G. Bryant and W. van Megen, J. Chem. Phys., 2007, 127, 
084505. 

54 N. Koumakis, M. Laurati, S. U. Egelhaaf, J. F. Brady and G. Petekidis, 
Phys. Rev. Lett, 2012, 108, 098303. 

55 P. A. Smith, G. Petekidis, S. U. Egelhaaf and W. C. K. Poon, Phys. Rev. 
E, 2007, 76, 041402. 

56 J. C. Conrad, H. M. Wyss, V. Trappe, S. Manley, K. Miyazaki, L. J. Kauf- 
man, A. B. Schofield, D. R. Reichman and D. A. Weitz, /. Rheol., 2010, 
54, 421^438. 

57 H. Konig, R. Hund, K. Zahn and G. Maret, Eur. Phys. J. E, 2005, 18, 
287-293. 

58 T. Hamanaka and A. Onuki, Phys. Rev. E, 2007, 75, 041503. 

59 P. Yunker, Z. Zhang and A. G. Yodh, Phys. Rev. Lett, 2010, 104, 015701. 

60 F. Ebert, P. Keim and G. Maret, Eur. Phys. J. E, 2008, 26, 161. 

61 L. J. Bonales, F. Martinez-Pedrero, M. A. Rubio, R. G. Rubio and F. Or- 
tega, Langmuir, 2012, 28, 16555. 

62 I. Biazzo, F. Caltagirone, G. Parisi and F. Zamponi, Phys. Rev. Lett, 2009, 
102, 195701. 

63 R. S. Farr and R. D. Groot, J. Chem. Phys., 2009, 131, 244104. 

64 N. Koumakis, A. B. Schofield and G. Petekidis, Soft Matter, 2008, 4, 
2008-2018. 

65 C. Mayer, E. Zaccarelli, E. Stiakakis, C. N. Likos, F. Sciortino, A. Mu- 
nam, M. Gauthier, N. Hadjichristidis, H. Iatrou, P. Tartaglia, H. Lowen 
and D. Vlassopoulos, Nat. Mater., 2008, 7, 780. 



10| 



